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Abstract Presently seven nearby radio-quiet isolated 
neutron stars discovered in ROSAT data and character- 
ized by thermal X-ray spectra are known. They exhibit 
very similar properties and despite intensive searches 
their number remained constant since 2001 which led to 
their name "The Magnificent Seven". Five of the stars 
exhibit pulsations in their X-ray flux with periods in the 
range of 3.4 s to 11.4 s. XMM-Newton observations re- 
vealed broad absorption lines in the X-ray spectra which 
are interpreted as cyclotron resonance absorption lines 
by protons or heavy ions and / or atomic transitions 
shifted to X-ray energies by strong magnetic fields of 
the order of 10^^ G. New XMM-Newton observations 
indicate more complex X-ray spectra with multiple ab- 
sorption lines. Pulse-phase spectroscopy of the best stud- 
ied pulsars RX J0720.4-3125 and RBS 1223 reveals vari- 
ations in derived emission temperature and absorption 
line depth with pulse phase. Moreover, RX J0720.4-3125 
shows long-term spectral changes which are interpreted 
as due to free precession of the neutron star. Modeling 
of the pulse profiles of RX J0720.4-3125 and RBS 1223 
provides information about the surface temperature dis- 
tribution of the neutron stars indicating hot polar caps 
which have different temperatures, different sizes and are 
probably not located in antipodal positions. 

Keywords stars: neutron • stars: magnetic fields • 
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1 Introduction 

After the discovery of RXJ1856. 4-3754 as an isolated 
neutron star in ROSAT data Walter et al..l99&.Walter fc 
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Il997j) . six further objects with very similar properties 
were found. Their X-ray emission is characterized by 
a soft, blackbody-like continuum which is little atten- 
uated by photo-electric absorption by the interstellar 
medium, indicatin g small dis tances (of the order of a 
few 100 DC. lPosselt et al.l2006j) . For the brightest object, 
RXJ1856. 4-3754, t his was soon confirmed by the par- 
allax measurement ll Waited 1200 it iKaolan et alJ 120023 
IWalter fc Lattimerll2002tl . Including new HST observa- 
tions the most recent revision of the parallax yield s a 
distance of 16lt\l pc l|van Kerkwiik fc KaDla'nll2006j) . A 
preliminary value for the parallax of RX J0720. 4-3125 is 
presented byivan Kerkwiik & Kaplan (2006} which cor- 
responds to a distance of 330^80° PC- 

No indication for a hard, non-thermal X-ray emission 
component was found in the X-ray spectra. Compared to 
normal radio pulsars also no strong radio emission is de- 
tected: While for RX J0720.4-3125 and RX J0806.4-4123 
deep radio observations failed to detect a radio counter- 
part (Joh nstonl2003tlKaplan et" .l2003bD . for RBS 1223 
and RBS 1774 the detection of we ak radio emission was 
claimed (jMalofeev et al.lf2006blal) . The constant X-ray 
flux on time scales of years, no obvious association with 
a supernova remnant and relatively high proper motion 
measurements for the three brightest objects suggest that 
we deal with a group of cooling neutron stars with ages 
of around 10^ years. Recent reviews which summarize 
the observed properties of the seve n thermally e mitting 
INSs k i iown toda y can be found in ^Treves et alJ l|200gi . 
iMotchI 1I2OOID and lHaberll (12004120051) . 

The discovery of the seven neutron stars (which are 
often called the Magnificent Seven, hereafter M7) with 
purely thermal X-ray spectra raised wide interest by the- 
oreticians and observers as promising objects to learn 
BfeffitiKnitifrs yspheres and the internal structure of neutron 
stars fe.g. lPaerelslll997j) . The apparent absence of non- 
thermal processes which would hamper the analysis of 
the X-ray emission allows a direct view onto the stellar 
surface. Limits on the radius of the star RXJ1856.4- 
3754 were used to co nstrain the equa t ion of state of 
neutron star matter l|Pons et all l2002t iTriimper et"al] 
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Table 1 X-ray and optical properties of the Magnificent Seven 



Object 


kT 


Period 


AmpUtude 


Optical 


PM 


Ref. 




cV 


s 


% 


mag 


mas/year 




RXJ0420.0-5022 


44 


3.45 


13 


B = 26.6 




1 


RXJ0720.4-3125 


85-95 


8.39 


8-15 


B = 26.6 


97 


2,3,4,5,6 


RXJ0806.4-4123 


96 


11.37 


6 


B > 24 




7,1 


RBS 1223("' 


86 


10.31 


18 


nisoccd = 28.6 




8,9,10,11 


RXJ1605.3+3249 


96 


6.88? 


? 


B = 27.2 


145 


12,13,14,15 


RXJ1856.5-3754 


62 




<1.3 


B = 25.2 


332 


16,17,18,19 


RBS 1774*'') 


102 


9.44 


4 


B > 26 




20,21,22 



= IRXS J130848.6+212708 

= IRXS J214303. 74-065419 
References: 

(1) Haberl et al. f2004a) (2) Haberl et al. (1997) (3) Cropp er et al. (20011 

(4) Haberl et al. (2004b) (5) dc Vrics ct al. (2004) (6) Motc h et al. (20031 

(7) Haberl & Zavlin (2002) (8) Schwopc ct al. (1999) (9) Hambarvan et al. (20021 

(10) Kaplan ct al. (2002a) (11) Haberl ct al. (2003) (12) Motch ct al. (1999) 

(13) Kaplan ct al. (2003a) (14) van Kcrkwiik ct al. (2004) (15) Motch ct al. (2005) 

( 16) Walter fc Matthews ( 1997') (17) Walter fc Lattimer (2002) (181 Burwitz et al. (2003) 

(19) van Kerkwiik fc Kulkarni (2001b'' l (201 Zampieri eTaTTMoi') (21) Zane erar(2005T (22) Komarova (20061 



I2004j) . The ROSAT PSPC spectra with low energy res- 
olution were consistent with Planckian energy distribu- 
tions with blackbody temperatures kT in the range 40 
- 100 eV and little attenuation by interstellar absorp- 
tion. More recent high resolution ob servations of the two 
brightest objects RX J1856.4-3754 (Burwitz ct al."200l") 
and RX J0720.4-3125 (.P aerels et al. .2001: Pavlov et al. 
l2002HKaDlan et alJl2003b^ were performed using the low 
energy transmission grating (LETG) aboard Chandra 
and the reflection grating spectrometers (RGS) of XMM- 



sorption. Apart from this absorption line with an equiv- 
alent width of 160 eV the X-ra y spe ctrum of RBS 1223 is 
blackbody-like. IHaberl et ~al\ l)2003|) also suggested that 
changes in the soft part of the X-ray spectrum of the 
pulsar RXJ07 2 0. 4-3 125 with pulse phase reported by 
CroDDcr et al. (200 j) are caused by variable cyclotron 
absorption (,Haberl et alJ I2004b1 . Cyclotron resonance 
absorption features in the 0.1-1.0 keV band are expected 
in spectra from strongly magnetized neutron stars with 
field strengths in the range of 10^° - 10" G or 2 x 10" - 



Newton. In particular the high signal to noise LETG 2 x 10^* G if caused by electrons or protons, resp ectively 
spectrum from a 500 ks Chandra observation of RX J1856.4~{see e.g. lZane et al ] l200lHZavlin fc Pavlov^l2002^ . Varia- 
3754 did not rev eal any significa,nt deyiation from a black- 
body spectrum l|Burwitz et alJl200lL 12003^. The statis- 
tical quality and energy band coverage of the RGS and 
LETG spectra of RX J0720.4-3125 are however insufH- 
cient to detect subtle narrow features in the spectrum 



tion of the magnetic field strength over the neutron sta r 
surface leads to a broadening of the line l)Ho & La^2004^ . 



ijPaerels et al.ll200lH . 

Four of the M7 exhibit clear pulsations in their X- 
ray flux, indicating the spin period of the neutron star 



The relatively long spin periods in the range of 3.45 s 
— 11.37 s (the majority of radio pulsars has periods less 
than 1 s) were the first indicator for strong magnetic 
fields. If the stars were born with millisecond spin pe- 



(RXJ0420.0-5022,RX J0720.4-3125,RXJ0806.4-4123and].iods^age estimates from neutron star cooling curves 



RBS 122 3) and there is ev idence that also RBS 1774 is 
a pulsar ijZane et alJl2005l) . A possible candidate period 
for RXJ1605. 3-1-3249 needs future confirmation. A sum- 
mary with blackbody temperatures derived from (phase- 
averaged) EPIC-pn spectra, spin periods, amplitudes for 
the flux modulation, optical brightness and measured 
proper motions is given in Table ^ together with key 
references for eac h object. Despite extensive searches in 
ROSAT data (e.g. iRutledge et alJ2003HChieregato et alJ 
l2005t lAgiieros et al.ll200(it) no further neutron star with 
similar characteristic properties was found. 

The first significant deviations from a bla ckbody spec- 
trum were reported bv lHaberl et al. I (|2003*) from EPIC 
spectra of RBS 1223. A broad absorption feature near 
300 eV was interpreted as proton cyclotron resonance ab- 



l)Page e t al. 200^ and proper motio ns (tracing the neu- 
tron star back to a likely birth place, iMotch et al.l l2006') 
of typically 10^ years require B fields of the order of lO^'^ 
G to decelerate the rotation of the stars to their current 
periods. This was recently confirmed by the accurate de- 
termination of the pulse period derivative P of 0.698 x 
10-" s s-i and 1.120 x 10" " s s"^ for RX J0720.4-3125 
and RBS 1223, respectively iKaplan fc van Kerkwiikl2005alb^ . 
In the magnetic dipole braking model this yields charac- 
teristic ages of 1.9 x 10^ years and 1.5 x 10^ years and B 
field strengths of 2.4 x 10" G and 3.4 x 10" G, respec- 
tively. Such strong fields supports the picture in which 
the broad absorption lines - if due to cyclotron resonance 
- at least for RX J0720.4-3125 and RBS 1223 are caused 
by protons. 
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2 Pulse-phase averaged X-ray spectra 

Significant deviations from the absorbed blackbody spec- 
trum were reporte d from XMM-Newt on EPIC observa- 
tions of RBS 1223 jHaberl et alJl200.^) and RXJ0720.4- 
3125 ijHaberl et al.l2004b^ and from XMM-Ncwton RGS 
spectra of RX J1605.3-h3249 ( van Kerkw ijk ct al. 2004). 
Also non-magnetic neutron star atmosphere models (e.g. 
iGansicke et al ■l2002HZavlin & Pavlovl200l fail to repro- 
duce the spectra. Iron and solar mixture atmospheres 
cause too many absorption features in particular at en- 
ergies between 0.5 and 1.0 keV which are not seen in 
the measured spectra. On the other hand the spectrum 
of a pure hydrogen model is similar in shape to that 
of a blackbody but results in a much lower effective 
temperature which would predi ct a far too high opti- 
cal flux fsee iPavlov et al.lll99(1^ . The modehng of the 
X-ray spectra was significantly improved by including 
a broad absorption line with Gaussian shape. Line en- 
ergies around 300 eV were found for RX J0720.4-3125 
and RBS 1223 from the medium energy resolution EPIC 
spectra, while a line with higher energy of 450 - 480 eV 
was discovered in the high energy resolution RGS spec- 
trum of RXJ1605. 3-1-3249. 



Spectral fits to the EPIC-pn spectra of the third 
brightest M7 star RX J1605. 3-^3249 obtained in 2003, 
January 17 and February 26 (satellite revolutions 569 
and 589) also considerably improve when adding a Gaus- 
sian absorption line to the absorbed blackbody model 
(see Table |21 and Fig. 0). However, the quahty of the 
fit with a reduced of 2.39 (aU XMM-Newton spec- 
tra presented here were obtained using the analysis soft- 
ware SAS version 6.5) is still not acceptable and worse 
than typical values around 1.5 one finds for joint fits to 
the available EPIC-pn spectra of RX J1856.4-3754 (see 
Fig. EJ or RXJ0720.4-3125. Adding more lines further 
improves the fit quality and with three absorption lines 
an acceptable is found. The best fit parameters for 
the models with different number of absorption lines are 
listed in Table El The width of the lines was fit as a 
single common parameter in the case of multiple lines 
(ct = 87 eV for the model with three lines). A remark- 
able result is the ratio of the line energies for the model 
with three lines. With E2/E1 = 1.46±0.02 and E3/E1 
= 1.94±0.06 the line energy ratios are consistent with 
1:1.5:2. It should also be noted that the depth of the 
lines decreases with line energy by a factor of about 5 
(in terms of absorbed line photon flux) from one line to 
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Table 2 Model fits to the EPIC-pn spectra of RX J1605. 3+3249 



Nh 


kT 


El 


Ea 


E3 


Fi 


F2 F3 


EWi 


EW2 


EW3 


Xred 




eV 


eV 


eV 


eV 


10-* 


photons cm-'^ s^^ 


eV 


eV 


eV 




0.15 


97 


















4.38 


1.1 


92 


445 






-67 










2.39 


1.4 


92 


440 


1.5Ei 




-75 


-43 








1.75 


2.04±0.04 


91.6±0.3 


403±2 


589±4 


780±24 


-43±1 


-8.0±0.8 -1.6±0.4 


96 


76 


67 


1.39 




0.3 0.5 1 

Channel Energy (keV) 



Fig. 2 EPIC-pn spectra of RX J1856.4-3754 obtained from 
five observations between 2002 and 2006 with the same in- 
strumental setup (small window readout mode and thin fil- 
ter). A joint fit with an absorbed blackbody model yields a 
reduced of 1.55 and demonstrates the stability of detector 
and source to better than 1%. The dip in the residuals around 
0.4 keV is a remaining calibration problem. 




0,5 

Channel Energy (keV) 



Fig. 3 RGS spectra of RX J1605. 3-1-3249 from the two ob- 
servations in 2003, January 17 and February 26 (satellite 
revolutions 569 and 589). The histograms represent the best 
fit EPIC-pn model with three absorption lines allowing only 
a re-normalization to account for possible cross-calibration 
problems. 



the next. A fit with the model with three absorption lines 
to the RGS spectra (only allowing a re-normalization fac- 
tor between the instruments) shows that also the RGS 
spectra are consistent with that model (Fig. O . 

Similar to RX J1605. 3-1-3249, new XMM-Newton ob- 
servations also indicate a more complex X-ray spectrum 
for RBS 1223. With a total exposure time of more than 
100 ks the increased statistical quality of the EPIC-pn 
spe ctra requires two abs orption lines for an acceptable 
fit llSchwope et alJ 120061) . With 230 eV and 460 eV the 
line energies are a factor of two apart. Also, in the case 
of RXJ0806. 4-4123 a simple absorbed blackbod y model 
yields an unacceptable fit to the EPIC-pn spectra l|Haberl 
l2004aD . A model with one line results in xled ~ ^ "^^ 
which is formally acceptable. Adding another line with 
an energy twice that of the first one, further improves 
the fit to xled — 1-05, but more higher quality data is 
required to confirm the significance of this result. 

Table 13 summarizes the magnetic field estimates for 
the M7 neutron stars utilizing either the magnetic dipole 
braking model for the pulsars with measured pulse pe- 
riod derivative P (B = 3.2 x 10^^(P x Py/^) or assuming 
that the absorption line is due to cyclotron resonance by 
protons (B = 1.6 x 10^^E{eV)/{l - 2GM/c^Rf/^). If 
multiple absorption lines were included in the spectral 



modeling, the lowest line energy was used for the esti- 
mate of the magnetic field. 

The brightest M7 star, RX J1856.4-3754, does not 

show pulsations with a very low ainplitu de limit of 1.3% 

I Ransom et al.l2002HBurwitz et al ]200l. However , Ivan Kerkwiik fc KuU 
( 2001al) discovered the existence of a Ha nebula around 
the star with cometary-like morphology aligned with the 
direction of its proper motion. Assuming that the neb- 
ula is powered by magnetic dipole braking and an age 
of the neutron star of ~5 x 10^ years which is inferred 
from its proper motion and the distance to the likely 
birth place , a m agnetic field st r ength of ^1 x 10^'^ G 
l erived llKapl an ct al.1 l2002d iBraie fc Romanil l2002t 



irumper et al.l I2004D . Such a field strength is consis- 
tent with the non-detection of a proton cyclotron feature 
which would have an energy below the sensitive range of 
X-ray instruments. 

In two cases (RX J0720.4-3125 and RBS 1223) accu- 
rate P measurements exist. For both stars also absorp- 
tion lines were found in their X-ray spectra which allows 
a comparison of Bjb (field estimate from dipole brak- 
ing) and Bcyc (field estimate from proton cyclotron reso- 
nance). The ratio Bcyc/Bdb is 2.33 for RX J0720.4-3125 
and 1.35 for RBS 1223 (for the spectral model with two 
lines). A similar low ratio for RXJ0720. 4-3125 would be 
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Table 3 Magnetic field estimates 



Object 


dP/dt 




Bdb 


Bcyc 




10-13 ss-i 


eV 


10^3 G 


10^3 G 


RXJ0420.0-5022 


<92 


330 


<18 


6.6 


RXJ0720.4-3125 


0.698(2) 


280 


2.4 


5.6 


RXJ0806.4-4123 


<18 


430/306('') 


<14 


8.6/6.1 


RBS 1223 


1.120(3) 


300/230'"' 


3.4 


6.0/4.6 


RXJ1605.3+3249 




450/400('') 




9/8 


RXJ1856.5-3754 










RBS 1774 


<60('') 


750 


<24'''' 


15 



Spectral fit with single / two lines 

With single hne / three hnes at 400 eV, 600 eV and 800 eV 
''^^ Estimate from Ha nebula assuming that it is powered by magnetic 

dipole breaking (KaDlan^^^^002c; ^r^^^^Romani^OO^; ^umge^^^^004) 
Radio detection: Malofeev et al. f 2006b ) 



obtained if a second (lower- energy) line around 140 eV 
exists. Unfortunately, the detection of a line at such low 
energies is outside the capability of the XMM-Newton 
instruments. 



3 Neutron star surface temperature 
distributions 

The first XMM-Newton observations of RX J0720.4-3125 
revealed spectral changes with pulse phase expr essed as 
hardness ratio variations ijCropper et al.ll200H) . Using 
polar cap models first constraints on the polar cap sizes 
and viewing geometries could be derived. The pulse pro- 
file in the 0.12 keV to 1.2 keV energy band was found 
to be approximately sinusoidal with a peak-to-peak am- 
plitude of about 15% with the hardness ratio being soft- 
est slightly before flux maximum. Spectra from differ- 
ent pulse phases indicated a temperature variation and 
a change in the low-energy attenuation, which - when 
modeled with a simple absorbed blackbody spectrum - 
appears as var iable ab sorption column density. Includ- 
ing more data iHaberl et al.l l|2004bh interpreted this as 
broad absorption line wit h variable depth. Pul se phase 
spectroscopy of RBS 1223 l)Schwope et al.ll2005[l also re- 
vealed that changes in temperature and absorption line 
depth are responsible for the spectral variations with 
pulse phase. The double- humped X-ray light curves of 
RBS 1223 in different energy bands were modeled by 
ISchwopc et al. {2QQ5) assuming Planckian radiation from 
a neutron star surface with inhomogeneous temperature 
distribution. Again, a simple model with two hot spots 
with temperatures Tf = 92 eV and T^ = 84 eV and 
full angular sizes of 8° and 10° which are separated by 
~ 160° can reproduce the data. A more physical tem- 
pe rature distribut i on ba sed on the crustal field models 
by iGeppert et alJ l)2004j) which can produce relatively 
strong temperature gradients from the magnetic poles 
to the equator was equally successful. 



The neutron star RXJ0720. 4-3125 is unique among 
the M7 as it shows a gradual change of its X-ray spec- 
trum over years which is accompa nied by an energy- 
dependent change in the pulse profile l|de Vries et al .12004 
IVink et alj r2004). Figure 21 shows folded light curves to- 
gether with hardness ratios for six different epochs. They 
were derived from observations with identical instrumen- 
tal setup to avoid systematic differences between differ- 
ent filters and CCD readout modes. As was found by 
Ide Vries et all (|2004j) the pulse profile became deeper 
with time with a pulsed fraction of ~8% in the year 2000, 
increasing to ~15% by the end of 2003. Moreover, the 
spectra became considerable harder during these years 
as can be seen in the increase of the average hardness 
ratio in Fig. H iHaberl et all l|200fi,) found first evidence 
for a trend reversal of the spectral evolution and inter- 
pret the involved period of 7.1±0.5 years as precession 
period of the neutron star. 

The overall hardening of the spectra until 2004 and 
the trend-reversal can also be seen from the pulse-phase 
resolved spectra which are shown in Fig.[Slfor four differ- 
ent XMM-Newton observations. Both, pulse-phase vari- 
ations and long-term spectral changes can be modeled 
by variations in the blackbody temp erature and absorp- 
tion line depth ijHaberl et alJl2006[l . In Fig. the hne 
equivalent width is plotted versus temperature kT. The 
temperature variation ('^2.5 eV) was smaller during the 
first observations and increased to ~6 eV, almost as large 
as the long-term change of ~8 eV seen in the phase- 
averaged spectra. This supports the idea that the tem- 
perature changes are caused by geometrical effects and 
different areas on the neutron star surface come into our 
view. In contrast, the amplitude in the equivalent width 
variation is ^40 eV and did not change much over the 
years. The relatively large variation in the line depth 
with pulse phase and the strong increase over the years 
indicates that the hotter polar cap which came into (bet- 
ter) view - possibly due to precession - plays the major 
role in the line absorption. More detailed analyses are 
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Fig. 4 Folded EPIC-pn light curves of RX J0720.4-3125 in two diflterent energy bands (soft S: 120 eV - 400 eV; hard H: 
400 eV - 1 keV) together with the hardness ratio. Different observations, all performed with the same instrumental setup 
(full-frame mode with thin filter), are shown from top left to bottom right: 2000 May 13, 2002 Nov. 6, 2004 May 22, 2005 
April 28, 2005 Sep. 23 and 2005 Nov. 12. For a detailed summary of all XMM-Newton observation in the years 2000 to 2005 
see Table 1 in Haberl et al. (2006). Pulse phases were calculated using the X-ray timing ephemerides ("All Data" solution) 
from lKaplan fc van Kerkwiik (.2005a) . 
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0.2 0.5 1 0.2 0.5 1 

Channel energy (keV) Channel energy (keV) 



Fig. 5 Pulse-pliase resolved EPIC-pn spectra of RX J0720. 4-3125 from four different observations. During the pulse the 
spectra vary as is shown by the different colours used for different phases (0.0-0.2: black, 0.2-0.4: red, 0.4-0.6: green, 0.6-0.8: 
blue, 0.8-1.0: light blue; with phase 0.0 defined as intensity maximum as in Fig. 0J. The intensity scale is linear to better 
resolve the changes at low energies. 



required to see if this is mainly a temperature effect (for ROSAT observations (which are subject to cycle count 
higher temperatures atomic line transitions are expected ambiguity) compared to the orig inal values published in 
to be reduced because of a higher ionization degree) or iKaplan fc van Kerkwiild (|2QQ^J. The reason for the dif- 



a vievsfing effect due to the dependence on the angle be- 
tween the direction of radiation propagation and mag- 
netic field lines. 

Since the spectral obs ervations do not ye t cover a 
complete precession period lHaberl et all l|20^fi^ also per- 
formed a preliminary analysis of published and new pulse 
timing residuals (from three new XMM-Newton obser- 
vations). This includes archival ROSAT, Chandra and 
XMM-Newton observat ions which are distributed over a 
total of ^12 years (se e 'Kapl an et al.ll2 002b^. 'Zane et al.' 
l2002HCropper et al.l2 004: Ka plan fc van Kcrkwiik 2005a) 
A sinusoidal fit to the phase residuals yields a period of 
7.7±0.6 years (Fig.Cj) consistent with th at derived from 



ferent periods might be due to the energy dependence 
of the pulse profile which can lead to systematic dif- 
ferences between instruments with different spectral re- 
sponse and/or the deviations from a sinusoidal shape of 
the light curves (which moreover evolve with time ) as it 
was ass umed in the analysis of IKaplan fc van Kerkwiild 
l|2005a|) . Further monitoring of RX J0720.4-3125 and a 
more sophisticated analysis and modeling of the data 
are required to confirm the precession model and further 
constrain its parameters. This can provide independent 
information about the interior of the neutron star and 
its distortion from a spherical shape e = (I3 — Ii)/Ii = 



4 X 10 whic h is larger than that re- 



the spectral analysis. In a similar a.na.lvsis lWKerkwiik fc i<Pittm jorrachopulsar (e.g. [ Jones fc Andersson| | 200n 
ll200a infer a shorter period of 4.3 years. Apart from us- Akgun et al.| |2006[) but smaller than that for Her X-1 
ing new Chandra data instead of the new XMM-observation A^^'^^^P^'^ al. || 1986|; [ Ketsaris et alj | 2000„) . 
they also used different phase residuals for the early 
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13-05-2000 (rev 0078) 
06-11-2002 (rev 0533/534) 

22- 05-2004 (rev 0815) 
28-04-2005 (rev 0986) 

23- 09-2005 (rev 1060) 
12-11-2005 (rev 1086) 




kT (eV) 

Fig. 6 Equivalent width of the absorption line vs. temper- 
ature kT derived from the EPIC-pn full-frame mode obser- 
vations with thin filter. Diamonds denote the values derived 
from the phase-averaged spectra. During the pulse the param- 
eters evolve counter-clockwise, the circle marks phase 0.0-0.2. 
This fi gure is a colour representation of Fig. 4 in .Habcrl ot aL 

ii2ro. 
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5x10* 5.1x10* 5.2x10* 
MJD (days) 



5.3x10* 



Fig. 7 Phase residuals for RX J0720. 4-3125. The black data 
points are reproduced from Kaulan & van Korkwijk (2005a) 
while red points are from the analysis of iHaberl et al.l (20061), 
which includes a re-analysis of published EPIC-pn results for 
a consistency check as well as three new data sets (the last 
three points). 



1 J ■ and long; periods |Cainilo et alJl200{]l: iMorris et al.ll2002l: 

^^^'^1 12006|) usedapola^p model similar to iMcLaughlin et a l.''2003') show! that radio emission can 

that applied bv lSchwope et all |2(}0ai to RBS 1223. In 



eluding free precession of the neutron star to explain the 
long-term spectral variations of RX J0720. 4-3125 they 
can at least qualitatively reproduce the variations in the 
X-ray spectra, changes in the pulsed fraction, shape of 
the light curve and the phase-lag between soft and hard 
energy bands. Like for RBS 1223, where the intensity 
peaks in the pulse profile are not separated by 0.5 in 
spin phase, the spots on the surface of RX J0720. 4-3125 
are probably also not located exactly in antipodal po- 
sitions. This can explain the observed evolution of the 
light curve and hardness ratio with precession phase. 



4 Discussion 

The new results from optical and X-ray observations of 
the Magnificent Seven, thermally emitting isolated neu- 
tron stars, strongly support the model of cooling neutron 
stars with ages around a million years. Large proper mo- 
tions measured for the three brightest objects make ac- 
cretion from the interstellar medium to re-heat the neu- 
tron star very inefficient. Absorption features in the X- 
ray spectra of most of the Magnificent Seven and first 
measurements of the spin period derivative from the two 
pulsars RX J0720.4-3125 and RBS 1223 consistently point 
to magnetic field strengths of the order of 10'^'^ G to 10^^ 
G. This places them at the long spin period and high 
magnetic field end of the radio pulsar distribution. How- 
ever, the P measurements still indicate field strengths 
below those of the magnetars (Fig. ^ . The discovery of 
a few radio pulsars with similar magnetic field strengths 



still occur at inferred field strengths close to the "quan- 
tum critical field" B^,. = mlc^/eh ~ 4.4 x lO^^ q There- 
fore it remains unclear if the M7 exhibit no radio emis- 
sion at all or if we do not detect them because their radio 
beam is very narrow due to their large light cylinder ra- 
dius and therefore does not cross the Earth. 

In a ~10^'^ G magnetic field cyclotron resonance lines 
at soft X-ray energies are expected to be caused by pro- 
tons or highly ionized atoms of heavy elements. The ra- 
tios of line streng ths in consecutive h armonic lines scales 
with Ecyc/(TOC^) ("Pavlov et al.'l980) and for the involved 
high particle masses no harmonic lines should be de- 
tectable. There is now more and more evidence for addi- 
tional absorption lines in the X-ray spectra of the Mag- 
nificent Seven. An explanation for the additional (or even 
all) lines can be atomic bound-bound or bound-free tran- 
sitions. In high B fields atomic orbitals are distorted into 
a cylindrical shape and the electron energy levels are 
similar to Landau states, with binding energies of atoms 
strongly increased. For hydrogen in a magnetic field of 
the order of 10^^ G the strongest atomic transition is 
expected at energy E/ eV ^ 75(l-H0. 131n(Bia))+63Bi3, 
with Bi3 = B/10" G l|Zavlin fc Pavlo v 2002). For the 
line energies found in the spectra of thermal isolated neu- 
tron stars this requires similar field strengths to those de- 
rived assuming cyclotron absorption. Atomic line transi- 
tions are expected to be less prominent at hi gher temper- 
ature s because of a higher ionization degree ijZavlin fc Pavlov! 
l2002|) . For a more detailed discussion of expected line fea- 
tures that can be pro duced in a hydrogen atmosphere see 
ivan Kerkwiik fc Kapl an ( 2006.) . It is also not clear how 
important the contribution of heavier elements in the 
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Fig. 8 P — P diagram for radio pulsars (black dots), mag- 
netars (blue dots) and the Magnificent Seven. Red dots show 
the direct P and P measurements for RX J0720.4-3125 and 
RBS 1223, while the open diamonds mark the P values ex- 
pected from the magnetic field estimates assuming proton 
resonance as origin of the absorption lines. For RXJ0720.4- 
3125 and RBS 1223 both field estimates are connected with 
a vertical line. Marked in green are radio pulsars which are 
also detected at high energies (incomplete for the millisecond 
pulsars) . 

atmosphere is b ecause of the strong g ravitational strat- 
ification forces l)Mori fc HailevI 12006^ . Further, the re- 
markable harmonic-like energy spacing of multiple lines 
needs to be understood. Studying their pulse phase de- 
pendence which can probably best be done for RBS 1223 
might shed light on this. 

The modeling of pulse profiles in different X-ray en- 
ergy bands is a powerful tool to unveil t he surfa ce tem- 
perature distribution of the neutron stars l)Zane fc TuroUa 
120061) . First studies of RX J0720.4-3125 and RBS 1223 
yielded constraints on the polar cap geometry. It is re- 
markable that for both stars asymmetric configurations 
of the hot polar caps were inferred. In both cases dif- 
ferent temperatures and sizes of the hot spots are found 
with the smaller spot being the hotter one. Probably, 
also in both stars the spots arc not located in exactly 
antipodal positions. This could point to magnetic field 
configurations of an off-centered dipole or the involve- 
ment of higher-order multi-pole components. Studies of 
the effects of a strong magnetic field on the temperature 
distribution in the neutron star crust show that configu- 
rations which include dipolar poloidal and toroidal com- 
ponents can i ndeed reproduce the observ ed temperature 
distributions (iGeppert et alJl2004 12001 . 
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